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Abstract. Measurements of charmed particle cross sections at HERA in the 
photoproduction and deep inelastic regimes are reviewed. The status of the comparison 
with perturbative QCD calculations is discussed. 

1. Introduction 

Open charm production at HERA has developed as a precision testing ground for fixed 
order (FO) NLO pQCD calculations as well as for other QCD inspired models. The 
bulk of the available data is on inclusive charmed particle production cross sections 
for several species. I will focus on reviewing the status of the comparison of the data 
with FO NLO pQCD calculations trying to identify the open issues for HERA Run II. 
The description of the kinematics for these observables requires two variables for the 
charmed particle momentum (an integration on the azimuthal angle of the charmed 
particle is generally assumed) and one (W) or two (Q^ and y, for instance) additional 
variables needed for the description of the inclusive photoproduction (PHP) or deep 
inelastic (DIS) cross section, respectively. The pseudorapidity, r] = —log(tan{9/2)), 
where 9 is the polar angle with respect to the proton (p) direction, and the transverse 
momentum with respect to the beam line, p^, are typically chosen as the charmed 
particle variables. Open charm production is intimately related to the gluon parton 
density in the p because the Photon-Gluon Fusion (PGF) process, sketched in Figure |1] 
(a), is expected to be one of the dominant production mechanisms. Since the invariant 
mass of the partonic final state, M|^^^, is significant due to the presence of two heavy 
quarks (M^^^^ > 4m^), the kinematics are such that the fraction of the p momentum 
carried by the gluon, Xg = XBj{Ml^j.^ + Q^)/Q^ rather than equal to xbj, as in the quark 
parton model (QPM) type of process. The sensitivity to the gluon density, fg, and to 
as can be read from the factorization formula for the PGF process written in equation 
(1), where a is the calculable partonic cross section. 

a^^XBj, Q^) oc e^asifir) ^fgixg, fif)a{xg, fif, fir) ; x™" = Xj; /"'^t^' (1) 

Higher order corrections and the fact that the cross section is measured in restricted 
regions of the charmed particle phase space complicate this simple picture. 

HERA Run 1 is already finished but the full data samples are just beginning to be 
exploited. HI and ZEUS have collected > 100 pb^^. This luminosity allows for the 
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selection of large charm samples. This is illustrated in Figure (c), from |T|, where 
~ 27000 D* mesons were reconstructed via the widely used K27r, D*^ —>■ (D^ 
K~7r+)7r+(+ c.c), decay channel 0, ||, ^ ||, |^. This decay mode has a low 
combinatorial background and a sizable branching ratio (0.0262±0.0010). D* mesons 
have been tagged also via the KAtt, D*^ — »• {D^ — »• i^'~7r"'')7r+7r"'"7r~(+ c.c), mode [|, 0] 
which has a factor of two larger branching ratio (0.051±0.003) but a larger combinatorial 
background. There are also data on Dg and mesons [|], Another window to 
study open charm is the identification of e~ coming from semileptonic decays [|l^ . This 
approach profits from the large branching fraction /(c e^) = 0.095 ± 0.009. 
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Figure 1. (a) Sketch of the PGF process and (b) charm excitation diagram, (c) D* 
peak from . 



2. Open charm photoproduction 

Real photon (7)-^ interactions have a hadronic component in which the 7 fluctuates 
into a hadronic final state (resolved). Thus, factorization has to be applied also for the 
7 and pdf's have to be measured. In addition to the direct PGF component, a large 
resolved contribution comes from the " charm excitation" process. Figure |l| (b). 

The ZEUS Collaboration § has measured the x^^^ = EjetsEre-"^ /{2yEe) 
distribution using D* dijets. The separation between resolved and direct is unique 
at LO, becoming scheme dependent at NLO. The comparison of the ZEUS data with 
the HERWIG LO predictions. Figure || (a), yields ~ 40% resolved component in the 
kinematic region of the measurement, dominated by charm excitation. Figure |^ (b) 
displays the comparison of the dijet data with a massive NLO pQCD prediction produced 
by the FMNR program |jl2[, which underestimates the low x^^^ values (resolved). In 
the massive approach only light quarks and the gluon pdf's are present in the p and the 
7. Another approach is to neglect the mass of the charm quark [1^ and treat it as a 
massless quark. The resummation of logs in Pt/^c results in a charm pdf. The main 
questions in open charm photoproduction are what is the relative weight of the available 
production mechanisms and what is the more appropriate model in the HERA range, 
massive or massless. However two caveats limit the precision of the test: It is known that 
the 7 pdf's, obtained from fits to Fg^, are poorly known in the HERA region and a new 
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iteration including PHP HERA data would be very welcomed. Moreover, pQCD cannot 
predict the fragmentation from a charm quark into a hadron. Thus hadron fractions 
and fragmentation functions are taken from e^e~ processes. The Peterson fragmentation 
function P{z) oc z~^[l — 1/z — e/(l — is convoluted with the charm quark 



distribution produced by the partonic matrix element, scaling the momentum of the 
charm quark to that of the D* meson by z. Note that this model is not Lorentz 
invariant because of the masses involved, /(c D*^) = 0.235 ± 0.007 ± 0.007 ||15| and 
are typically used for D* cross section predictions. 
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Figure 2. Untagged PHP D* dijet cross section differential in x^^^ compared to 
HERWIG (a) and FMNR (b). Inclusive D meson differential cross sections da/dp^ 
(c ) and da/df]^ (d) where D stands for Dg (dots) and D* (triangles). Thick (thin) 
curves in (b) (c) (d) are FMNR predictions with parameters set to rric = 1.5(1.2) GeV, 



p^{0.5y/ml +p^) and fip = 2(4)^^^. 



ZEUS has measured D* p and Ds P untagged PHP, where the scattered lepton 
escapes undetected through the beam hole, using 38 pb^^ of ^/s = 300 GeV data. 
The kinematic region of the untagged measurements is: 130 < W < 280 GeV, 
< 1 GeV^ (Qmedian ~ ^ ■ 10~^ GcV^ ). The Dg mesons are reconstructed using 
the — > 07r"*" — >■ {K~^K~)tt~^{+c.c.) channel. Experimental restrictions result in the 
phase space limitation 3 < pt{Ds) < 12 GeV and | rj^Dg) \< 1.5. 

The fraction of the Dg to the D* cross sections in the above phase space region 
is measured to be fj§f| = 0.41 ± 0. 07 {stat.)1:^-^i{syst.) ± 0.10(6r.). This quantity 
is sensitive within the string model to the strange suppression parameter, 7^. 7^ 
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= 0.27 ± 0.05 ± 0.07(6r.) is obtained with the help of JETSET, being in agreement 
with e^e~ and therefore supporting the universahty of c fragmentation found in the 
earher HI pubhcation [Q, and the more recent ZEUS results [IC, l7]. 

Figures (c) and (d) display differential cross sections of Ds and D* in untagged 
PHP. Most of the branching ratio error, which dominates the 7^ measurement, can 
be neglected in the comparison with the calculation since it scales the data and the 
prediction in the same way. The curves are the FO NLO pQCD prediction from 
FMNR using MRSG (p, a,=0.114) and GRV-G HO (7, a,=0.111) pdf's and the 
Peterson fragmentation parameter e=0.035. Normalization uncertainties coming from 
the luminosity measurement (1.7%) and the hadronization fractions (3% for D*and 9% 
for Dg) are not plotted. The Dg data support the same conclusions of the more accurate 
D* data, namely: 

• The normalization of the theory is too low. 

• The shape of the r] distribution is not well described. 

Massless predictions give a better description of the data both in shape and 
normalization as can be seen in Figure |^. The parameters for the calculation are pdf 



CTEQ4M (p, as = 0.116), e = 0.1, = ^ml+pl ; m, = 1.5 GeV ; iip = '^f^R- The 
difference between the curves is the 7 pdf. Surprisingly enough, the predictions using 
GS-G HO, with a c pdf equal to the u pdf, are the closest to the data. The better 
agreement of the massless calculations is not expected since the resummation effects are 
small for pt(-D*)<10 GeV, i.e. in the region of the measurement 
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ZEUS has also measured tagged PHP ^ yielding a similar level of description by 
FMNR. These data are better described by the BKL model [0, which fits the untagged 
sample. 

The HI Collaboration has measured open charm PHP selected using the 44 (33) 
m taggers 0. The kinematic region of the measurement is: < 0.009(0. )GeV^, 
<W >= 88(194) GeV, Pt{D*) > 2(2.5) GeV and | rapidity(L'*) |< 1.5. The data are 
compared in Figure § with the FMNR prediction using pdf's MRSTl (p, 0^=0.1175), 
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GRV-HO (7), e = 0.035, ^ir = Jml+pl 



1.3 - 1.7 GeV ; /if = 2/1^. This 



calculations used an old value of /(c D*), 15% higher than the latest result (0.235). 
Also the use of a larger value of as than in the case of the ZEUS analysis produces higher 
values of the prediction. The HI Collaboration finds good agreement with FMNR and 
uses this calculation to extract Xgg{xg). 
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Figure 4. Tagged PHP D* cross section differential in rapidity (a) (b) and px (c). 



3. Open charm electroproduction 

Earlier HERA {y/s = 300 GeV) and fixed target data (^i ~ 30 GeV) show that open 
charm electroproduction is dominated by PGF, thus providing direct sensitivity to the 
gluon pdf in the p. FO NLO pQCD calculations are available in the form of a MC 
integrator (HVQDIS) [|20|. Only light quarks and the gluon are present in the p as 
pdf 's. Since open charm electroproduction is sensitive to the gluon, this process allows 
a test of factorization and in particular of the universality of the gluon pdf, comparing 
the data with NLO predictions which use as an input the pdf's obtained from fits to 
inclusive F2. The same fragmentation limitations mentioned before for PHP apply here. 

ZEUS has measured D* production in DIS using £ ~ 37 pb~^ at i/i = 300 GeV in 
the kinematic region I < < 600 GeV^, 0.02 <y< 0.7, 1.5(^.5) < Pt{D*) < 15 GeV 
and I ri{D*) |<1.5 0. The differential cross sections are shown in Figure ^ The 
two decays used, D* — > K27t{K4tt), are in agreement as can be seen in the pt{D*) 
distribution. The settings of the NLO HVQD IS predict ion are: pdf ZEUS NLO fit 
(a,=0.119), e = 0.035 (in the lab), fi^ = fip = ^Am^^ + Q\ = 1.3 - 1.5 GeV. An old 
value of /(c D*^), 9 % lower than the latest result (0.235), was used. The comparison 
of the data with this prediction (open band) yields the following observations: 

• Good agreement with HVQDIS in Q^, xbj and W. 

• HVQDIS is shifted with respect to the data in the t]{D*), x{D*) distributions. 

However, the shaded band, obtained by folding JETSET and the same NLO calculation, 
describes better the data. The NLO calculation is folded with JETSET by means 
of reweighting the pt{c),ri{c) two dimensional distribution of a LO MC after the hard 
interaction and letting JETSET hadronize. A shift into the forward direction due to the 
string model is observed. A similar picture is obtained using HERWIG or ARIADNE. 
Figure ^ shows that the shift in rj coming from the full fragmentation model (PS+string 
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Figure 5. Differential cross sections for D* production from the K2tx final state 
(solid dots) and from the KA-k channel (open triangles). The open band corresponds 
to the standard Peterson fragmentation function. The shaded band shows the NLO 
reweighted JETSET MC. x{D*) (f) is defined as 2p*{D*)/W, where p* is in the 7*-p 
CMS frame. 



model) is similar to the effect coming from the string model alone. This shift is expected 
in models which include the interaction of the colour charges in the final state since the 
colour configuration is similar to e"^e" — qqg^ where the string effect was observed. In 
PGF processes the remnant is acting as a colour octet which drags the c, c into the 
forward region. 

ZEUS has measured cross sections of e~ coming from semileptonic decays of charm 
in DIS using £ ~ 34 pb^^ ||TI|. This measurement relies on the dE/dx deposition of 
the tracks in the CTD to separate the e~ signal from the hadronic background. The 
kinematic region is: 1.2 < Pe < 5 GeV, | r/e |< 1-1, 1 < < 1000 GeV^, 0.03 <y < 0.7. 
The method is complementary to the D* analysis since, on one hand, it profits from 
the larger branching fraction to go to higher and, on the other hand, it probes 
harder c quarks, providing a pt enriched sample. The price to pay is that it covers 
a smaller region of the c quark phase space. The data are in agreement with NLO 
HVQDIS, which in this case uses the settings: pdf GRV94 HO (a^ = 0.111), e = 0.035, 



l^R = l^F = Y^^^c + ^'^d rUc = 1.3 — 1.7 GeV. 

ZEUS has also measured D* production cross sections integrating two samples 
at different i/i, C ^=300,318 ~ 45,38 pb^^ The integrated cross section in the 
kinematic region 10 < < iqoo GeV^, 0.04 < y < 0.95, 1.5 < Pt{D*) < 15 GeV 
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Figure 7. Scmilcptonic e differential cross sections. 



and I r]{D*) |< 1.5 is a(e+p e+D*^X) = 2.33 ± 0.12(stat.)i[{;J^(syst.) nb. This 
result and the differential cross sections compare reasonably well with NLO HVQDIS 
using three different pdf's, see Figure ||. Although the measurement limits have been 
pushed towards high and xbj, MC studies show that there is no sensitivity to intrinsic 
charm (IC) yet. Another issue in this region is whether variable ffavour number schemes 
(VFNS), which include a c pdf arising from the resummation of large loglQ"^ /mc) can 
give a better description of the data than that of the FO NLO pQCD predictions. Little 
difference between FFNS and VFNS's was found in |^ until Xsj > 0.1 and high Q^. 
However, a slight difference appears at ~ 10 GeV^, Figure 0. Even if the effect is not 
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Figure 8. D* cross sections differential in logio{Q^) (a) and logio{xBj) (b). 



an artifact introduced by the use of HVQDIS to produce VFNS predictions differential 



in tlie c phase space or by the normahzation of all curves at 



the question 



is whether it can be distinguished from fragmentation and/or mass uncertainties which 
are of the same order. 



Figure 9. FO NLO pQCD 
(labeled Exact) is compared with 
two VFNS predictions and with 
HI and ZEUS data. Figure taken 
from I2I. 
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The HI Collaboration has performed a measurement of D* production in DIS 
using £ ~ 18 pb~^ 0]. The integrated cross section in the kinematic region 1 < 

Q 



< 100 GeV^, 0.05 



< y < 0.7, 1.5 GeV < Pt{D*) and | r]{D*) 
= 8.50 ± 0.42(stat.)lJ:?^(syst.) ± 0.65(mod.) nb. 



< 1.5 is: 

<^kin[g'p cu'^A.) = 8.0U ± u.42(stat.j_:_5;7g(syst.j ± u.DS^moa.j nD. This can 
be directly compared with the ZEUS measurement [0 in the similar kinematic region 



1 < < 



600 GeV^, 0.02 < y < 0.7, 1.5 < Pt{D*) < 15 GeV, | r]{D*) \< 1.5: 
(TKiN{e+P e+D*^X) = 8.31 ± 0.31(stat.)i[JJ(syst.) nb using HVQDIS to interpolate 
between the two measurements. 
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which is compatible with 1. This, and the comparison of the differential cross sections 
in logio^Q"^), see Figure]^, support the fact that HI and ZEUS data agree within errors. 

Hi's data are compared with HVQDIS using pdf GRV98 (a,=0.114), mc=1.3- 
1.5 GeV. In addition to the standard Peterson (in YP CMS) with e=0. 035-0.1, a 
transverse momentum with respect to the charm quark is given to the D* meson, 
according to the function exp{—apf) with < p^ >^350 MeV. This results in a prediction 
of 7.02 — 5.17 nb. The HI analysis also compares the data with CASCADE p2[, which 
implements the CCFM evolution. In this model an unintegrated gluon density, fitted to 
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Hi's F2, and mc=1.3-1.5 GeV is used. The fragmentation is performed by JETSET with 
6=0.035-0.1 yielding 10.77 — 8.04 nb. The measured single and double differential cross 
sections are compared to the same calculations in Figure |10[ CASCADE shows better 
agreement in normalization, while the previous HVQDIS prediction is too low at forward 
Tj, low zd* (Figures ^ (e) and (f)). The double differential cross section (Figure |10| (g)) 
shows that the disagreement with HVQDIS is concentrated at low p-i{D*) and forward 
ri{D*), which is correlated with low zd*. Moreover there is a poor agreement with 
CASCADE at high pTiD*). 
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3.1. Extraction of F^'^ 

F2'^ can be defined in terms of the double differential cross section integrating in the 



dxdQ'^ 



{[1 + (1 - ?/)2]Ff (x, g2) _ y2^f }. Several 



full c quark phase space: 
assumptions are made in the extraction: 

• The contribution of Ff^ (< 1%, according to the FO NLO calculations) is neglected. 

• Bound charm (<2.5-4.5 % [0) is neglected. 

• /(c ^ £)*+) from e+e^ is vahd at HERA. 

• The extrapolation factors, f^j^? are well described by the model used to 
extrapolate. 

The relevance of the last assumption is illustrated by the size of the extrapolation factors, 
which go from <2 at high to <4 at low for D* data and can be as large as <20 at 
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low for e data. The actually applied procedure to extrapolate outside the kinematic 
region (KIN) in px?''? is: 



The extrapolation is not strongly affected by the fragmentation models because the 
procedure is not very sensitive to the shape of the r] distribution. The information 
contained in the Fl'^ plots in Figure |ll] can be classified according to its dependence on 
the extrapolation: 

(i) Since F2^rneasl ■^2^theo ~ '^kW/'^kin coustructiou (see eq. |^), the comparison 
of -Fl'^meas "w-ith the -Fl'^^/ieo used to extrapolate represents the level of 
agreement at the cross section level and is independent of the extrapolation. The 
semileptonic double differential cross section a']^J]^{Q'^ ,XBj) is well described by 
HVQDIS (GRV94) as can be seen in Figure [l^ (a). The same is true for the D* 
data and HVQDIS (ZEUSNLO), Figure [Tl| (b). On the other hand, good agreement 
is observed between the unextrapolated HI data and CASCADE, Figure |lT] (c). 

(ii) Statements about the shape of -F2'^^eas model dependent: Fg"^ rises steeply as 
XBj decreases and, when plotted vs. at constant xsj 0, does not scale in the 
measured region. This can be interpreted as coming from the rise of the gluon via 
the PGF dominance of open charm electroproduction. ¥2'^ represents ~ 25—30% 
of F2 for Q'^> 11 GeV^ and low XBj- Figure |Tl| (b) has two regimes: it flattens at 
low XBji where both ¥2^ and F2 are dominated by PGF and, neglecting all mass 
effects, ¥2'^ IF2 ~ 2_i_ 2_i_''2_i_ 2 1 i ~ 4/11 ~ 0.36. The denominator of the ratio, F2, 
rises at high xbj due to the valence contribution, forcing the ratio to drop. 

(iii) Comparing the shape of F^^^^^s "with other models, or two F^^^^^s extracted 
assuming different theories is, to a large extent, a comparison between two models, 
rather than with data. Of this kind are all comparisons that can be made in 
Figure |TI| (d): F!^^^^^^ from three independent samples extracted with HVQDIS 
(GRV 98, GRV 94, ZEUS NLO) seem to be compatible. The prediction using the 
HI fit agrees reasonably well with the shape of HVQDIS (GRV 98, GRV 94, ZEUS 
NLO). 

4. Summary 

The overall picture is that of reasonable agreement among all data samples, obscured by 
the different interpretation of these data, especially concerning the level of description 
by FO NLO pQCD calculations. 

In PHP, HI observes an adequate description by FMNR and uses these data to 
extract Xgg{xg). However, ZEUS observes a deficit in the FMNR prediction, specially 
in the forward region. Alternative models/explanations are massless calculations, BKL, 
fragmentation effects or large NNLO corrections. 

The agreement within theoretical uncertainties of the FO NLO pQCD HVQDIS 
predictions using pdf extracted from DGLAP fits to the inclusive F2 scaling violations 
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with the HI and ZEUS D* production cross section and the ZEUS semileptonic cross 
section provides support for the universahty of the parton distributions of the p and the 
vahdity of factorization in charm electroproduction. However, HI finds better agreement 
with CASCADE predictions than with those of HVQDIS. 

5. Outlook 

Measurements of open charm production will benefit from the expected increase of a 
factor of 5 in luminosity, together with the instrumental improvements of the HI and 
ZEUS detectors in HERA Run II. This can result in a extended coverage of the phase 
space, especially in the forward region, which could be used to: 

• Improve the understanding of charm photoproduction. 

• Increment the sensitivity to IC in DIS, since forward rj correspond to higher XBj- 

• Look for indications to distinguish between the VFNS approach and FO NLO. 

• Minimize the extrapolation dependence of Fg'^. In this case it would be useful to 
explore channels beyond the usual D* modes, whose lower px limit is fixed because 
of the good correlation between pt{D*) and ptItTs). 

Ongoing studies on c-meson spectroscopy fl], |10|, [l^, diffraction ||23| or double tagged 
final states will benefit from the increase in statistics. 

On the theoretical side, an improved treatment of fragmentation, matching the NLO 
calculation to MC fragmentation models as JETSET, could lead to a better description 
of the data. It would be also desirable to identify a more distinctive signature for CCFM 
evolution. The final goal is to include open charm cross section data in global pdf fits. 
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Figure 11. Measured Fj^^ from semileptonic data ||Tl| vs. xbj at fixed values (a) 
compared to the model used for the extrapolation (HVQDIS, GRV94). Ratio of the 
measured Ff from D* data extrapolated with (HVQDIS, ZEUSNLO) over F2 vs. 
XBj at fixed values (b). Measured from D* data [Q vs. XBj at fixed values 
compared to the model used for the extrapolation (CASCADE) (c). The data of (a) 
and (b) compared with the HI data B extracted with HVQDIS,GRV98 (d). 
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